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In the framework of a field theoretical approach we study Maier-Saupe nematogenic fluid in contact with a hard 
wall. The pair interaction potential of the considered model consists of an isotropic and an anisotropic Yukawa 
terms. In the mean field approximation the contact theorem is proved. For the case of the nematic director 
O ' being oriented perpendicular to the wall, analytical expressions for the density and order parameter profiles 

^ ! are obtained. It is shown that in a certain thermodynamic region the nematic fluid near the interface can be 

+^ i more diluted and less orientationally ordered than in the bulk region. 
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Due to orientational ordering, nematic fluids near a surface show richer behavior than in the case of 
simple fluids. Among them are the anchoring phenomena, whereby the surface induces a specific orien- 
tation of the nematic director with respect to the surface (lfl. In order to understand this phenomenon, 
during the past decade the Henderson-Abraham-Barker (HAB) approach , previously developed in 
the theory of isotropic fluids in contact with solid surfaces Q|, has been employed. In this approach, the 
description of the fluid density profile reduces to the solution of the Ornstein-Zernike (OZ) integral equa- 
, tion for the fluid particle-wall distribution function calculated from the known fluid particle distribution 

function in the bulk. In the framework of the HAB approach, the application to the bulk of the nematic 
model, analytically solvable at the level of the mean spherical approximation (MSA) |5Q, makes it possible 
to investigate the role of orientational-dependent molecular interactions with the surface in anchoring 
phenomena QSl. However, in the MSA, this approach does not take into account the contribution from 
long-range molecular interactions and, as a result, does not satisfy the exact relation known as the con- 
tact theorem . According to this theorem, the contact value of the particle density near a hard wall 



for a neutral fluid is determined by the pressure of the fluid in the bulk volume. 



Recently, the density field theory, previously developed for ionic fluids near a hard wall , has 
been applied to the description of simple fluids with Yukawa-type interactions near a hard wall Hill . 
In both cases, the developed approach yielded correct results. In this theory, the contributions from the 
mean field and from the fluctuations are separated. In II ill , it was shown that the mean field treatment 
of a Yukawa fluid near a wall reduces to solving a non-linear differential equation for the density profile 
while the treatment of fluctuations reduces to the OZ equation with the Riemann boundary condition. 

In this paper, the density field theory, developed in for simple fluids at a hard wall, will be gen- 
eralized to nematic fluids at a hard wall. To this end, we consider Maier-Saupe nematogenic fluid model 
(T3>[T3] as one of the simplest models that account for the isotropic-nematic phase transition. For simpli- 
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fication we consider a fluid of point uniaxial nematogens interacting through the pair potential 

v(ri2,QiQ 2 ) = v (ri 2 ) + v 2 (ri 2 )P2(cos0i2), (1) 

where the first term vo {r\2l - [Aq I r^) exp (- ocq Tyi) describes the isotropic repulsion and the second term 
with v 2 (ri 2 ) = [A2lr\2) exp (-0:2 12) describes the anisotropic attraction between particles (Aq > 0, A2 < 
0), r\2 denotes the distance between particles 1 and 2, Q. — [8, <p) are orientations of particles, P2 (cos 612) = 
(3 cos 2 0i2 - l)/2 is the second order Legendre polynomial of the relative orientation 612. 

The application of the density field theory to the description of bulk properties of such nematic 
fluids was considered in (3. El- It was shown that beyond the mean field approximation, the repul- 
sive isotropic term in Q} is very important for the description of the nematic phase. Within the field- 
theoretical formalism, the Hamiltonian is a functional of the density field and can be written as the sum 
of entropic and interaction terms 

BH[p{r,Q.)] = J p(r,n)[ln(p(r,Q)A fl A|,)-l]drdfi (2) 

+ ^ J v(ri2,QiQ2)p(ri,ni)p(r2,Q2)dridr 2 dQidQ2, 

where B - II k^T is the inverse temperature, dQ = {l/4n) smOdOdcp is the normalized angle element, 
p(r,fi) is particle density per angle, so that / p(r,Q)dQ = p(r), Aj is the thermal de Broglie wavelength 
of the molecules, the quantity A^ 1 is the rotational partition function for a single molecule 11611 . 

In this paper, we will restrict our consideration to the mean field approximation (MFA) which is the 
lowest order approximation for the partition function. In the canonical formalism it corresponds to fixing 
the Lagrange parameter A, so that the following relation is true for the singlet distribution function 



8p(r,Q) 

As a result, 



(3) 



where 



p(ri,ni) = p bulk (ni)exp|-ig f v(n 2 ,Oi0 2 ) [p(r 2 ,Q2)-p bulk (fi2)]dr 2 dn 2 J, (4) 

1 

p b ^(n) = p b exp[-(a:|Sb/af)P 2 (cos0)]/ J dcos0 exp [- (afSb/af) P 2 (cos0)] (5) 



is the singlet distribution function for the bulk nematic in the MFA defined in the framework of 
the Maier-Saupe theory Qjl x\ - Anp\ s pA2, Pb is the bulk value of the fluid density, S b = 
(1/pb) fo P2(cos0)p bulk (Q)dcos0 is the bulk value of the orientational order parameter. 

In order to integrate with respect to the angle in (D, we should separate the angle Q ;n between the 
particles and the director and the angle Q, nw between the nematic director and the surface. To this end, 
we express the Legendre polynomial in the potential v(ri2,ni02) in terms of spherical harmonics Y2 m (fi) 
as P 2 (cosni 2 ) = (1/5)1 Y* m {a ln )Y 2m {a 2 n). As a result, 



p(Ti,n.i n ,n wn ) 



- expi- [ V lti,a w ^- v£] - WOi re ) \v 2m (.Ti,n wn ) - V^jl, (6) 

P U'ln) I v5 m I 

where the mean field potentials 

V o 0ti,a wn ) = P Jv (ri2)p(r2,Q wn )dr 2 , (7) 
V 2m {ri,n wn ) = B / v 2 (ri2)S2m(r2,O u ,„)dr 2 . (8) 
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The bulk values of these quantities are V b = x^la^, V 2Q = x^S^la 2 ,, V 2 h m = for m + 0, where x^ = 
4np h pA , 



p(r,Q. wn ) = j pOr,Qi n ,a wn )dSli n 

is the density profile. The property 

S 2m {r,n wn ) = -j= j p(r,D.i n ,D. wn )Y 2m {n ln )dD. ln = p(r,n wn ) sZ m (i,n. wn )> 



(9) 



(10) 



where S2 m (r,Q. wn ) are the order parameter profiles. Far from the wall S^fr, — ► Sb, S 2m {r,D. wn ) — ► 
for m + Q. 

Taking the gradient of equation (6) we have 

1 1 „ 

— — - — - — -Vp(r,Qi„,Q, tm ) = E (r,Q w „) + — 2^i2m(nin)E 2 m(r,n wn ), 

p\X,ll\ n ,ll wn ) ^5 m 

where we define an equivalent of the electric field by 

E (r,n wn ) = -VV (r,n wn ), E 2m (r,n wn ) = -VV 2m (r,n u ,„). 
Due to the properties of the Yukawa potential 

[A-a 2 ) V {r,n wn ) = -4n/3A p{r,D. wn ), 



(11) 



(12) 



(A-a|) V 2m (r, = -4ti p A 2 S 2m {r, Q wn ). 



(13) 
(14) 



Due to the translational invariance parallel to the wall, the functions considered depend only on the 
distance z to the wall, and replacing Si~3\ and dl4l into fPT) , we obtain 



_d_ 

dz 



p(z,n wn ) a 2 ) _J_ F 2 f n ) 

Pb 2x1 



E 



"0 

2 ^2mi z '^wn) ~ E 2m {z,Q, wn ) 
ylX 2 



1 



2% 



= 0. 



(15) 



In the bulk, when z — «■ oo, we have p{z,VL wn ) — ► pb, -Eo(-z,fi«,n) — * 0, E 2m {z,D. wn ) — ► 0, VbCZjOum) — * 
Vq , %(z,Qim) -► V^g and V 2m {z,D. wn ) —■ for m ^ 0. From Equation (IS) we see that the quantity in 
brackets is constant regardless of the angle D M between the director and the wall, and, therefore, it can 
be evaluated, for instance, in the bulk as the reduced pressure /3PI pb within MFA 11511 : 



BP xi x\ 

Pb 2a\ 2a\ b 

Outside the system, where there are no particles, we have another invariant which is simply 



Q v 2 (z,n wn ) - ^-E 2 (z,n wn ) + £ 



2x 



2x 



0*2 2 ^"2 

~S V 2m ( z '®-wn) ~ — 2 E 2m {z,D, wn ) 



(16) 



(17) 



its value is zero far from the interface and, therefore, it is zero at the interface as well. From the continuity 
of the potential and its derivative due to equation fP3l and (14), we see that this is also true at the wall 
just inside the system z = 0+. Thus, 



p{0 + ,n wn ) a% 2 



1 



Pb 



+ -^v Ho + ,n lvn )-- 1 E^o + ,n wn ) 



2x 



E 



2#n 



~2 ^2m^+'^i"i) Z~2^2m^+'^^n) 



p{0+,Cl wn ) 



(18) 
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Since this quantity is constant, we obtain the so-called contact theorem 

pp = P (o+,n wn ). 



(19) 



We should note that the contact theorem was usually proved for isotropic fluids near a hard wall 

The result obtained here is probably the first verification of the contact theorem for anisotropic fluids at 

a hard wall. 

From equations fT0)-fT4t we have a set of six differential equations for unknown functions 
p{r,D.i n ,D. wn ), S 2m (r,n wn ), E (r,n wn ), E 2m {r,Q. wn ), V {r,D. wn ), V 2m {r,n wn ). We note that in the case 
when the director is oriented perpendicular to the wall, D. wn = 0, the singlet distribution function is axi- 
ally symmetric. Consequently, in the equations considered only the terms with m = will be present. In 
this paper, we will restrict our further consideration to this special case and consider the solution of the 
obtained differential equations in the linear approximation. Far from the wall p(r,Q) — ► p bulk (Q). After 
linearization of expression l[6) 



p'(z,Q) = [Eolz) + E 20 {z)P 2 (cos6)] p buIk (n) 



(20) 



and we have the following system of equations 



p'(z) = [E (z)+S b E 2 o(z)]pb, 



S' (z) = 



Mz)S b + E 20 (z)-(Y£ ) a 
b 



Pb, 



V ' (z) = - E (z) , V 2 ' (z) = - E 20 (z) , 
E' (z) = - a 2 Q V (z) + {x\ I p b ) p (z) , 
E 20 {z) = -a 2 2 V 2a {z) + {x\l p b ) S 20 {z) , 



(21) 

(22) 

(23) 
(24) 
(25) 



where the prime denotes a derivative by z and (Y 2 k ) n = (1/pb) / F 2 q(Q)p (Q)dcos0 (fc = 1,2). This 

o 

system reduces to two second order differential equations for Eq[z) and E 2 q{z) 

4'(z) = Bb(z) («g + a 2 ) + E 20 (z)x 2 S b , 



E 2 ' (z) = E {z)x 2 S b + E 20 {z) [a 2 + {x 2 /5) (Y 2 ) n ) , 
which can be solved with the boundary condition that should include the contact theorem dT9l . Thus, 



(26) 
(27) 



p(z) 
Pb 



S 2 o(z) 

Pb s b 



= 1 



H{iY 2 Q )n- (Y 20 ) 2 n ) n ^_ AqZ 



x 2 S b 



Bie 



X 2_ a 2_l x 2 [{Y 2 Q)n _ {Y2a) 2 a) 
27 B2e 

x\ S b 



X 2 ^ Sb 



"2 



(28) 
(29) 



where 



l 0,2 



x n + a n +- 



(Y 2 l Q )n + a 2 2 ± 



2 2 
X + a o " 



(30) 



Bi 



Xy S b 



2[A 2 -X 2 ) 



i 2 X 2 - 

_Q_ _2_ 

y2 + 



(^/5)(F 2 2 > n 



B 2 



^2 $b 
2^|" 



(31) 



The values of parameters Aq and X 2 coincide with similar parameters obtained in the bulk case after 
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P(z)/p b 



Sb 


= 0.63 (T* 


= 2.41) 


Sb 


= 0.56 (T* 


= 2.56) 


Sb 


= 0.47 (T* 


= 2.67) 



S* 20 (zVS, 



zct„ 




Sb = 


0.63 (T* 


= 2.41) 


Sb - 


0.56 (T* 


= 2.56) 


Sb = 


0.47 (T* 


= 2.67) 



Figure 1. (Color online) The density and the order parameter profiles in the linearized approximation 
for Pb/ag ~ °- 5, A a IA 2 = 2.2, a la2 - J 
reduced temperature T* = l/Cp^o^)!- 



for Pb/«o ~ 0-5, A0/A2 = 2.2, aala2 = 1.25. Different lines correspond to different values of [i.e., the 



including the Gaussian fluctuations 11511 . Similarly to the bulk case, parameters Ao and A 2 characterize the 
screening of the repulsive isotropic and the attractive anisotropic interactions, respectively. The density 
profile p(z) and the order parameter profile S^iz) = S2o{z)/p{z) calculated from equations ll28l-(l29) are 
presented in figure [1] As we can see, both profiles have layer-like forms. At the surface, the fluid is more 
dense and orientationally ordered than in the bulk. But at higher distances there is a large region where 
the fluid is more diluted and less orientationally ordered than in the bulk. This result is obtained in the 
MFA. For isotropic fluids, as was shown in Hill , the inclusion of fluctuation terms leads to the depletion 
effect. We can suppose that a similar effect will take place for the considered anisotropic fluid. It means 
that the inclusion of fluctuation terms does not change our conclusion drawn from figure [1] concerning 
the existence of a more diluted and less ordered region of the fluid near the surface compared to the bulk 
region. 

Other interesting phenomena that can be observed within the framework of the formulated MFA 
approach are connected with the possibility that the angle between the nematic director and the surface 
can change near the surface. Consequently, phases with the order parameters S^iz) = Szmiz)! p{z) {m + 
0) can appear near the surface. These phenomena will be investigated in a separate paper. 

In conclusion, using the field theoretical approach we have formulated the mean field approximation 
as the first starting point for the description of a nematic fluid at a hard wall. For the first time, an exact 
derivation of the contact theorem for anisotropic fluids has been presented within the MFA. It has been 
shown that the contact value of the density profile is determined by the pressure of the fluid in the bulk 
volume and does not depend on the angle between the nematic director and the surface. For the case of 
the director being oriented perpendicular to the wall, in the linear approximation for the MFA, we have 
obtained analytical expressions for the density and order parameter profiles. It has been shown that at 
some values of the parameters of the model considered, the fluid near the surface can be more diluted 
and less orientationally ordered than in the bulk. 
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HeMaTMMHMM n/ii/iH 6i/i?i TBepfloT noBepxini y Ha6/ii/i>i<eHHi 
cepeflHboro no/in 

M. ro^oBKcP, I. KpaBLjiEP, fl,. fli KanpicP 

iHCTHTyT (|)i3HKW KOHfleHCOBaHHX CHCTeM HAH YKpaiHH, By/1. I. CBEHL^bKOrO, 1, 7901 1 /lbBiB, YKpaiHa 

J1a6opa"ropm e/ieicrpoxiMiT, xiMii noBepxoHb i eHepreTUHHoro MO,qe.moBaHHH : Biflfli/iem-iH xiMii Bumoi 
Hai;iOHa/ibHOi luko/ih napiTex, a6. CKpuHbKa 39, n/i. >Kyccio, 4, 75005 napnx, ct>pai-mjfl 

B paMKax TeopeTMKO-no/ibOBOro ri\p,xop,y BHB4aeTbCH HeMaToreHHnw n/inH Mafiepa-3ayne 6i/ia TBepflOi no- 
BepxHi. y po3MHflyBaHiii MOfle/ii napi-iufi noiei-mja/i B3aeMOflii cK/iaflaeTbca 3 i30TponHoro Ta aHi30TponHoro 
KDKaBiBCbKnx ,qo,qaHKiB. y Ha6/in>KeHHi cepeflHboro no/in flOBefleHO KOHTaKTHO TeopeMy. fl,nn BwnaflKy, ko/ih 
flnpeKTop HanpaB/ieHufi nepneHfli/iKy/iapi-io flo CTiHKH, OTpuMaHO aHa/ii™4Hi Bnpa3n fl/ifl npo<|>i/iiB rycTHHH Ta 
napaMeTpa nopaflKy. noKa3aHO, mp y neBHM TepM0flHHaMi4Hifi o6/iacri HeMaTUHHnw n/iwH no6/iw3y noBepxHi 
MO>Ke 6yTn 6i/ibLU po3piflxeHWM i Memu opieHTauji'lHO BnopflflKOBam/iM, Hix b o6'eMHi£i o6/iacri. 

KniOMOBi c/ioBa: HeMaTmHUM nni/iH MafiepaSayne, TeopemKO-nojibOBi/id niflxifl, noBepxHfi po3p,tny, 
KOHTaKTHa reopeMa 
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